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out of equilibrium with a time-dependent Hamiltonian H(t). A measure of nonequi- 
librium temperature for the open system evolving from an equilibrium is proposed 

as the time-averaged energy expectation value T(t) = Ti((H(t))^ / {I(t))^), where 
I{t), the action operator, satisfies the quantum Liouville-von Neumann equation and 
determines the true density operator. It recovers the result for the adiabatic (quasi- 
equilibrium) and nonadiabatic (nonequilibrium) evolution from one static Hamiltonian 
to another and takes into account the particle production due to the intermediate pro- 
cesses. 
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1 Introduction 



A system becomes nonequilibrium (out of equilibrium) when its coupling constants change rapidly 
through the interaction with an environment. In inflationary scenarios, the inflation of the Uni- 
verse would rapidly cool the inflaton and matter fields and this supercooling process would 
have happened far from equilibrium (see ref. l|). Another phenomenon out of equilibrium 
is the recently introduced preheating mechanism, in which the parametric resonance of bosons 
and fermions coupled to the oscillating inflaton leads to catastrophic production of particles 
|2], |3|, |, ||, . In quark-gluon plasma the second order phase transition out of equilibrium is ex- 
pected to play an important role in quark-confinement ]7], ||, |§. Recently the nonequilibrium phase 
transitions have been actively studied to understand the dynamical process of domain growth and 
topological defects (see ref. |T 



The quantum statistics, however, has not attracted much attention for these nonequilibrium 
systems though it is important in understanding the multi-particle phenomenon. Further a mea- 
sure of temperature has not been assigned to such a nonequilibrium system because temperature is 
widely believed to make a physical sense only for an equilibrium system. But a quantum statistical 
theory has been proposed for nonequilibrium system based on the assumption that the microscopic 
world is correctly described by the Schrodinger equation and the quantum Liouville-von Neumann 
equation [jllj. For the nonequilibrium system the frequently-used operator pH{t) = e~ H ^^ kT /Z H 



cannot be used since it is not the true density operator. Instead the true density operator that 
satisfies the quantum Liouville-von Neumann is given by another operator pi(t) = e~ I ^ t ^ kT /Zj, 
from which quantum statistics such as the Bose-Einstein or the Fermi-Dirac distribution is de- 
rived. Hence it may be suggested that the nonequilibrium temperature should be defined out of 
the Liouville-von Neumann operators I(t) and the dynamical energy operator H(t). 

In this paper a quantitative measure of temperature is proposed for the nonequilibrium system. 
The analysis will be confined to the effective theory of an open boson or fermion system whose 
Hamiltonian H(t) explicitly depends on time. The spatial homogeneity of the system will be 
assumed and the interaction among constituent particles will be neglected. The time-dependent 
nature of the system comes mainly from the interaction with an environment. Further, to simplify 
the analysis, we focus on the open boson or fermion oscillator with time-dependent frequency 
u(t). The time-dependent boson oscillator has been extensively studied as an exactly solvable, 
nonstationary, quantum system [|TJ, [L3|, [TJJ. Also an open fermion system has been studied as 



the nonequilibrium system JT5|. For the time-dependent oscillator we derive the nonequilibrium 
temperature using the exact density operator pi(t) and the Hamiltonian operator H(t) as giving 
the dynamical energy. Many solvable models are investigated, which evolve from one static 
Hamiltonian to another and have the limit of the adiabatic (quasi-equilibrium) or the nonadiabatic 
(nonequilibrium) evolution. 

The organization of this paper is as follows. In Sec. II the relation between temperature 
and energy expectation value in the high temperature limit is discussed qualitatively and the 
temperature relation is derived approximately for the adiabatic and nonadiabatic evolution. In 
Sec. Ill, using the action operator and the exact density operator, the nonequilibrium temperature 
is defined by the ratio of the time-averaged energy expectation value to the expectation value of 
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the action operator. The nonequilibrium temperature is shown to be the quantum analog of the 
classical one from classical distribution function. In Sec. IV the nonequilibrium temperature is 
applied to several models of physical interest and extended to the time- dependent fermion system. 
In the nonadiabatic evolution from one static Hamiltonian to another, the final temperature 
is obtained in terms of the initial temperature, initial and final frequencies, and the particle 
production rate. 



2 Adiabatic vs. Nonadiabatic Evolution 



We start with an ensemble of boson oscillators with the frequency u in a thermal equilibrium 
at temperature T. The thermal state of the system is described by the density operator pu = 
e~ H ' kT / Zh and has the energy expectation value 

where % is the Planck constant and k the Boltzmann constant. In the high temperature limit 
kT ^> hoo, the thermal energy expectation value has the limiting value 

(H) PH = kT. (2) 

This may suggest that the temperature of the system should be derived from the energy expec- 
tation value at least in the high temperature limit. 

We now evolve the system from one static Hamiltonian Hi to another Hf by adiabatically 
changing the frequency from Ui to uuf. If the system was in a thermal equilibrium, it would 
maintain quasi-equilibrium at each moment from the initial equilibrium to the final one. Both 
the initial and final states are assumed to have high temperatures. Then the ratio of the initial 
and final temperatures is given by the ratio of the corresponding energy expectation values: 



^> = ( H f)p f 



(3) 



where pi and pf denote the initial and final density operators pu- If the density operator pu does 
not change significantly, eq. (||) can be approximated by 

T i (Hi)* <* " 

The particle production due to the change of frequency (parameter) is negligible during the 
adiabatic evolution. 

In the nonadiabatic evolution, the rapid change of frequency drives the system out of thermal 



equilibrium. Such a nonadiabatic process necessarily involves the particle production |IBL [T7], ITS 
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The creation and annihilation operators at the initial and final times are related with each other 
through the Bogoliubov transformation 



a(t f ) = aa(U) + /3a f (*,•), a)(t f ) = a*a\u) + /3*a(U), 
where the coefficients satisfy 



(5) 



a*a-P*/3 = l. (6) 

It is interesting to compare the expectation values of the final and initial Hamiltonians with 
respect to the number states, \n,ti), of initial time t%, whose ratio is given by 



It therefore follows that 



{Hi)m 
(Hf) Pi 



^) x (1 + 2(3* (3). 



UJi 



^ j x (1 + 2(3*13). 



UJ; 



(7) 



(8) 



The factor (3* (3 is the number of particles created per unit volume and depends on the intermediate 
nonadiabatic process during the evolution. As will be shown later, eq. (Q) is related to the 
temperature ratio. If a device measures the temperature of the system with time-averaged energy 
density, the final temperature would be given by 



T f = T t x 



^(1 + 2(3* (3) 



(9) 



3 Nonequilibrium Temperature 

To define a quantitative measure of temperature that recovers eqs. (|J) and @ in the limiting 
cases, we study an open boson system described by a time-dependent oscillator 

H(t) = ^ + ^ q \ (10) 

where the frequency us(t) explicitly depends on time. The quantum states of the system are 
governed by the time-dependent Schrodinger equation 

ih^^- = H(tMq,t). (11) 

The statistical properties of the system are determined by the density operator that satisfies the 
quantum Liouville-von Neumann equation 

in^ + [p(t),H(t)]=0. (12) 

The physical quantities which can be directly determined by eqs. (0) and (|T2|) are (H(t))y and 
(p(t))i&, from which all other statistical variables such as the temperature and free energy will be 



derived. Here \ty) denotes any exact quantum state of eq. (|TT|) . 
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To find the exact quantum state of eq. (O), we use the idea of invariant operators first intro- 
duced by Lewis and Riesenfeld [IS], who found a quadratic invariant of position and momentum. 
In refs. |2D], |2T], [TT], [I3J, a pairs of linear invariant operators satisfying eq. (|TJ) are introduced 



a f (t) 



VK 



[<p(t)p - (p(t)q] 



a(t) = —[y\t)i>-<p*{t)q\, 

where ip is a complex solution to the classical equation 

(p{t) +u 2 (t)(p(t) = 0. 

The Wronskian condition 

ip{t)<p*{t)-tp*(t)ip{t)=i 
guarantees the standard commutation relation 

\a{t),a ] {t)} = 1. 



(13) 

(14) 
(15) 

(16) 



Hence the time-dependent operators ^(t) and a(t) play exactly the same role as the standard 
creation and annihilation operators. The Fock space of exact quantum states consists of the 
number states of a){t)a{t) [fH|| 



1 



1 



1/4 / \ (2n+l)/2 



q 



cxp 



where H n is the Hermite polynomial. The density operator satisfying eq. (|T2|) is given by 

1 



(17) 



where 



Pi(t) 



I(t) = TiUi 



-I(t)/kTi 



;is) 



(19) 



Here the free parameters T, and u>i will be fixed as the temperature and frequency, respectively, 
if the system starts from an initial equilibrium. 

The action operator, gives an energy analog (/(t))* at each moment, where \& denotes 
the Fock state (|T7|), whereas (H(t))^ is the dynamical energy. We may write the expectation 
value of the density operator as 



-{I(t))*/kTi 



-(H{t))y/kT(t) 



Zj Zj 
where the time-dependent temperature is given by 



T(t) = Tj x 



V(/(0)J' 



(20) 
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The process of taking the expectation value with respect to the Fock basis can be understood 
from the correspondence principle between quantum and classical theory. The classical density 
distribution is given by 



where I(t) is the classical action 



Pi(t) 



I(t) = UOi 



-i{t)/m 



/Zi 



a*(t)a(t) + 



1 



(21) 



(22) 



with the classical invariants [21 



a*(t) = -i[<p(t)p-<p{t)q] 
a{t) = [<p\t)p-<p\t)q\. 



(23) 



The density distribution also can be written as 

1 

PiW = ^exp 



where 



T c {t) = Ti x 



H(t) 
kT c (t) 

(H(t)\ 



(24) 



The temperature (|2~0|) is the quantum analog of the classical result (|24"D, provided that the energy 
of H(t) and I(t) is evaluated appropriately in phase space. 



A few comments are in order. First, the temperature relation 
is replaced by pi(t). From the expectation values 



still holds even when 



(H(t)) Pl 
(/(*)>« 



-coth| 



huji \ 

2kTi) 
Tiuj, 



[ip*{t)v{t)+u 2 {t)ip*{t)ip{t)l 



^icoth!-- 
2 \ 2kT, 



(25) 



we find the temperature at any later time 



T(t) (H(t)) 



pi 



Ti 



pi 



(26) 



Second, if the time scale of oscillations determined by u>(t) is much smaller than the thermal 
relaxation time of measuring device, the time average should be taken over the period: 



V (i(t)) 9 J 



(27) 



Roughly speaking, the nonequilibrium temperature is determined by the time-averaged energy of 
the system at each moment. 
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4 Applications 



We now consider the adiabatic or nonadiabatic evolution of boson oscillator from one static 
frequency uj,i to another ujf. The solution to eq. (|l4j) in general has the asymptotic form 

p—iu>it 

= (28) 



^M^> (29) 

where a and j3 determine the Bogoliubov transformation (|5|) and satisfy the condition (|6]). The 
final state (JTT|) , into which eq. (fffi) is substituted, oscillates with the period n/wf. Then the 
temperature (|27D , obtained by taking the time average of eq. (|6|), is given by 



Tf = Ti x 



^-(l + 2/3*/3) 



(30) 



The temperature ( p0|) is the same as the temperature @ of the nonadiabatic evolution, which 
was obtained qualitatively. It also reduces to the adiabatic result when \/3\ <C 1, that is, the 
particle production is negligible. So the temperature (|27j) recovers the nonadiabatic evolution 
from one static Hamiltonian in equilibrium to another. 

We investigate the models with a parameter adjusting adiabaticity. The first model has the 
frequency which changes according to the modified Poschl- Teller type |22| 

<™-* + i£&m- (31) 

It has two asymptotic regions at t — > =Foo with the same frequency ujq. The strength of interaction 
is determined by A(A — 1) and the rate and duration by r. The factor for pair production is given 
by 

Fp= f*^_\\ (32) 

In the adiabatic limit of large r (r ^> 1), the particle production is exponentially suppressed by 
e -2iru) T an( j result (|3|) is recovered. One prominent consequence is the temperature enhance- 
ment in the nonadiabatic limit of small r (r < 1), which leads to f3*(3 ^> 1. Another model with 
an adiabatic parameter is given by the frequency 



= + wj) " 2 ~ <"/) tanh ( Z I • (33) 
It also has two asymptotic regions at t = =Foo with frequencies u>i and Uf, respectively. It is found 



-W/)tanh^. 



23 



smn(7rru;ij sinn(7rra;/j 
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In the adiabatic limit of large r, the particle production is again exponentially suppressed by 
e -7ir(wj+w/-|wj-a;/|)_ j g wor fa noting that in both models any nonadiabatic process with different 
t leads to different final temperatures. 

The next model is the preheating mechanism of light bosons or fermions coupled to an oscil- 
lating inflaton || ||], |J ||, § . To simplify the analysis, we consider the model frequency of a scalar 
field mode 

uj 2 {t) =ujI + 6(t)tof cos cut, (35) 

where the parametric coupling to the inflaton turns on at t = 0. There are unstable narrow bands, 
in which the solution to eq. (|T4|) grows exponentially due to the parametric resonance 

,2 



Then the temperature (P7D is given by 



% x 



CO 



exp 



CO 



-iut + — * 
2co 



co% 



-<! — + — + 1 ^ 

CJ 



MM* 



(36) 



(37) 



ujq uj uoquo \ 2lo / 

The thermal energy of exponentially growing temperature comes from the decaying inflaton. The 
temperature ( p7|) correctly describes the growing temperature in the preheating mechanism. 

Finally, we extend the above analysis to the nonequilibrium system of time-dependent fermion 
oscillator. Pairs of fermions are also produced when the frequency changes in time JP, 18 . For 
an ensemble of fermion (b, ¥) and its antipartcle (d, <$) with time-dependent parameters, pairs of 
linear invariant operators satisfying Eq. (^) are found in ref . [|I5| : 

b(t) = a b {t)b + (3 b {t)d\ P(t)=a* b (tp + PZ(t)d (38) 
d(t) = a d (t)d + p d (t)b\ S(t) =a* d (t)S + [3* d (t)b. (39) 

These operators satisfy the anticommutation relations 

{S(t),S+(t)} = l, {d(t)J\t)} = l, (40) 

iOl) or (1271) holds with 



and all the other anticommutation relations vanish. The general result 
the new action operator 

I(t) = huOip{t)b(t) - d(t)$(t)}. (41) 

So, when the fermion system evolves from one static Hamiltonian Hi with uoi in an equilibrium 
to another Hf with uof, we arrive at the final temperature 



T f = T l x 



LOj 



(42) 



The above temperature, which is valid for the adiabatic (j3 = 0) and nonadiabatic (/3 ^ 0) 
evolution, is the fermion analog of the boson result (|30"|). The coefficient relations due to the 
anticommutation relations 



a, b, 



(43) 



limit the fermion production due to the Pauli-blocking. 
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5 Conclusion and Discussion 



We have studied the nonequilibrium evolution of a time-dependent boson or fermion system and 
proposed the nonequilibrium temperature through eq. fl20D or (p7|) . The quantum statistics is 
defined quantitatively and rigorously for the nonequilibrium system under the assumption of the 
Schrodinger equation (quantum law) and the quantum Liouville-von Neumann equation. For the 
time-dependent boson or fermion oscillator these two equations can be solved simultaneously in 
terms of quantum invariant operators, which in turn can be found in terms of classical solutions of 
motion. The true density operator is constructed from a particular quantum invariant I(t), called 
the action operator. Using the well-known fact that in the high temperature limit the energy 
expectation value of a thermal ensemble of boson or fermion oscillator is proportional to the 
temperature, we define the nonequilibrium temperature fl2~0|) or (p^J) from the dynamical energy of 
the Hamiltonian opeartor H(t) and that of the action operator 1(f) of the true density operator. 
It recovers the temperature (||) for the nonadiabatic evolution from one static Hamiltonian to 
another, not necessarily in equilibrium. The nonequilibrium temperature is applied to many 
physical models in Sec. VI. 

Finally we discuss whether the nonequilibrium temperature ( p0|) or ( pTf) can be applied to a 
system whose final static Hamiltonian has a lower temperature. In the lower temperature limit 
the vacuum energy dominates over thermal energy and leads to the expectation value 



(6h„ = ^- (44) 



Therefore it seems that the qualitative derivation of temperature in Sec. II fails in this model. 
However, the direct analysis below shows that the temperature ( p0|) or fl27f ) may still be applied 
to this model. As a temperature lowering model, we consider a second order phase transition, 
which is described by the sign changing frequency 

uj 2 (t) = e(-t)uoj-e(t)uj). (45) 

Here 9 is the step function which takes 6(t) = 1 when t > and 6(t) = when t < 0. The solution 




to eq. (|lj) is given by 



<ff{t) = cosh(uj ft) — i— smh(uj ft) 

x'-^'f [ ' u f 

The expectation value (H(t))^ with respect to any number state (|TT| ) vanishes because the kinetic 
and potential energies contribute equally with opposite signs The final temperature Tf 

approaches zero when the instantaneous quench process continues for a sufficiently long time. In a 
finite quench the final temperature needs not to be exactly zero. This implies that the temperature 
proposed in this paper may have a wide range of validity including the lower temperature region. 

Not considered in this paper is the process of the nonequilibrium system evolving toward the 
final equilibrium state with the true density operator p# = e~ H fl kT f / 'Z^. The effective theory of 
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time-dependent boson or fermion oscillator cannot describe this process properly since the final 
states oscillate rapidly. To get the non-oscillating final Fock states, the final solution to eq. (14) 
should have the form fl25p , which in turn leads to the initial solution of the form ( |2"9"D . This implies 
that our effective theory of the open system cannot describe completely the nonequilibrium process 
from one equilibrium to another. Even with the more elaborated model that describes correctly 
the evolution from one equilibrium to another, the particle production due to the nonadiabatic 
process leads to the chemical potential. These points together with statistical relations will be 
addressed in a future publication. 
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